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We report measurements on the synchronization properties of organ pipes. First, 

we investigate influence of an external acoustical signal from a loudspeaker on the 

sound of an organ pipe. Second, the mutual influence of two pipes with different 

^ . pitch is analyzed. In analogy to the externally driven, or mutually coupled self- 

O I sustained oscillators, one observes a frequency locking, which can be explained by 

_0 ' synchronization theory. Further, we measure the dependence of the frequency of 

'c/2 ' 

J>^' the signals emitted by two mutually detuned pipes with varying distance between 

the pipes. The spectrum shows a broad "hump" structure, not found for coupled 
^ ' oscillators. This indicates a complex coupling of the two organ pipes leading to 

Tij- ■ nonlinear beat phenomena. 
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I. INTRODUCTION 



- ^ , Sound production in organ pipes is traditionally described as a generator-resonator cou- 

' pling. In the last decades, research has been concerned with the complex aeroacoustic 



processes which lead to a better understanding of the sound generation in a flue organ pipe. 
The process of sounding a fiue-type organ pipe employs an airstream directed at an edge, 
the labium of an organ pipe. An oscillating "air sheet" is used to describe the situation in 



which the oscil 



pipe sound 



m 



ations of the jet exiting from the flue are responsible for the creation of the 



Using the "air sheet" terminology, it is pointed out that the oscillation 



is controlled not by pressure, as in earlier investigations 



mm, 



but by the air flow 
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The situation becomes more involved if two flue organ pipes are close in sounding fre- 
quency jLnd in spatial distance. Then a synchronization of the pipes, a frequency locking, 
I5L llfiL Il7l | . This has a direct importance for the arrangement of organ pipes in a 



occurs 



common Orgelwerk. The effect has been known byorgan builders for a long time and taken 



into account intuitively in the design of organs 



18l |. Some measurements of the acoustic 



field or of its dependence on the mutual coup 
been reported first in the early 20th century 



ing or on the distance between the pipes have 
Iq . but no theoretical explanation was given. 



In this article, we report acoustic measurements and give a possible explanation by means 
of modern synchronization theory [19]. To a certain approximation, an organ pipe can be 
considered as a self-sustained oscillator, explained in detail below. In this context, the syn- 
chronization of a pipe by an external, acoustical force is important, such that we carried out 
additional measurements, where a pipe stands aside a loudspeaker whose frequency is tuned 
around the pipes pitch. As well in this case, the pipe is found to be synchronized. 

As a side result of these measurements, we can clarify a discussion about the nature of 
the strong amplitude decrease for two coupled pipes, already observed by Lord Rayleigh [l^ 
and late, by ve.y detailed n,ea„t= of Boua^e M. As a possible ...te.p.etatio,, the 
so-called oscillation death has been given in jl9|, which means that there is an oscillation 
breakdown at the the pipe mouth, and all the energy is dissipated. From our results, we rule 
out such a scenario and suggest an antiphase oscillation which yields destructive interference 
of the emitted acoustic waves. This holds for two neighboring pipes, as well as for the pipe 
standing aside a loudspeaker. 

This article is structured as follows: In section |n] we explain briefly the generation of 
sound in an organ pipe and review basic concepts from synchronization theory. In section 
mil we report on detailed measurements of the synchronization of a loudspeaker positioned 
directly on the side of a pipe and observe two adjacent pipes. The results are interpreted 
within the frame of synchronization theory. We give an explanation of why a model of two 
oscillators works out so well, as indicated in Further, the dependence of the frequency 
spectrum on the distance between two detuned pipes is investigated for a fixed amount of 
detuning. Finally, we conclude with section HVl 
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II. BASIC PRINCIPLES 

A. Sound generation in organ pipes 

Sound generation in organ pipes has been repeatedly investigated 0, 12, |2^. Here, the 
beauty of musical sound generation is paired with complex aerodynamical phenomena; their 
coupling to the acoustic field has been understood to a reasonable degree in the last 30 years 
(see the review 

Let us consider a single organ pipe: The wind system blows with constant pressure 
producing a jet exiting from the pipe flue. Typical Reynolds numbers, corresponding to a 
free jet are of the order of 10^, depending on the pressure supplied and the pipe dimensions, 



see 



12^ 



2l|. The jet exiting from the flue generates a pressure perturbation at the labium of 
the pipe which travels inside the pipe resonator and is reflected at the end of the resonator. 
This pressure wave returns after time T to the labium where it in turn triggers a change of 
the phase of the jet oscillations. After a few transients, a stable oscillation of an "air sheet" 
at the pipe mouth is established. This oscillation of the wind fleld couples to the acoustic 
modes and a sound wave is emitted. Thus, the system can be considered as a damped linear 
oscillator (the resonator) which controls the periodic energy supply of the air-jet by the air 



column oscillation. We want to adapt this idea of an organ pipe as a se^ 



f-sustained oscillator 



2J| , or synchronization. 



23| to understand the mechanism of mode- or frequency locking 
One feature of a self sustained oscillator is the occurence of oscillations even for constant 
driving. The physical mechanism is the balance of energy losses and energy input. The 
classical example is the van der Pol oscillator (Fig. QJ, where a Triode (or more modern, 
a tunnel diode) acts as "negative resistance" whose response is fed back to an electrical 
oscillator. In the case of an organ pipe, energy is supplied by the wind system, losses are 
due to internal dissipation in the spatially extended resonator, and radiation. The negative 
resistance is represented by the air stream at the pipe mouth and the feedback coupling 
is given by the boundary conditions between the "air sheet" wind fleld and the air-column 
oscillations in the resonator. It can be described by impedances j3| , detailed results for two 
coupled Helmholtz resonators are given in . Of course, the exact form of the limit cycle 
and thus the acoustic oscillations depend on the details of energy losses, as the quality factor 
of the resonators, the radiation of sound, and the understanding of the energy supply by 
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the coupling to the jet. The general mechanism, however, remains untouched, since only a 
nonlinear growth of losses and supply is needed for a self-sustained oscillator to work. Here, 
we do not want to investigate these details, but rather emphasize the general character of 
our investigations. 



feedback 




FIG. 1: Left:Schematic diagram for a self sustained oscillator. Right: resulting limit cycle behavior. 



The sound emitted by this complex system can be evaluated by the coupling of the 
flow and the acoustic modes in the framework of aeroacoustic modeling on the basis of 
Lighthills analogy. According to Howe the sound generation is dominated by 

the singularity at the edge of the labium. Numerical approaches still suffer from the very 
expensive runs needed to resolve the large range of excited scales and the proper choice of 
boundary conditions 



B. Synchronization of self-sustained oscillators 



In his book. Lord Rayleigh states "When two organ pipes of the same pitch stand side 
by side ... it may still go so far as to cause the pipes to speak in absolute unison in spite of 



inevitable small differences" jl^. He describes the so-called "Mitnahme-Effekt" (loosely the 
take-along effect), known by organ builders In Q the phenomenon is interpreted 

in the frame of synchronization theory. Hereafter we will use the terms frequency locking 
or synchronization to be synonymous with the "Mitnahme-Effekt". In the effect has 
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been analyzed heuristically, while Stanzial has investigated the dependence of the frequency 
locking on the detuning or frequency difference of two single pipes in by acoustic mea- 
surement. He modeled the effect by two coupled oscillators, without giving a physical reason 
of the coupling. A treatment of the mode-locking phenomenon, specific for musical instru- 
ments is very nicely given in [2^ . In the following we will rely on synchronization theory to 
show that any self-sustained coupled oscillator generically shows synchronization. 

Probably the most important feature of self-sustainment is the occurrence of an attracting 
limit cycle. It appears due to two properties: nonlinearity, and the energy balance between 
losses and driving. For a linear, damped oscillator a limit cycle solution does not exist, 
the only possible attractor is the trivial solution. Nonlinearity allows for the dependence 
of frequency on amplitude which constitutes a mechanism to drive the system towards an 
amplitude at which the regular oscillations are established. Since the amplitude corresponds 
directly to the mechanical energy of the system, this is right at the point of equality of losses 
and supply - the limit cycle. 

We want to sketch the principles for synchronization of an oscillator with an external 
driving and then explain the interaction of two oscillators. In the following we will discuss 
some equations in terms of the angular frequencies and return later to frequencies when the 
measurements are concerned. 

A Single Oscillator: Let us consider the equation for a harmonic oscillator with negative. 



nonlinear resistance: 



X = -UqX - ^f{x) = -ujIx - x^ , (1) 



where Uq is the angular frequency of the harmonic system and f{x) is a nonlinear function 
with at least one part of positive and another part of negative slope. Now we discuss 



the basic principles by the concrete example of 



;he van der Pol oscillator and give general 
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In the case of the van der Pol oscillator 



results without further derivation, for details see 
f{x) = ax — and = {LC)~^. The energy supply is accounted for by a, the losses 

by (3. If the inharmonicity is not too great, the amplitude on the limit cycle differs little 
from the harmonic oscillator and the phase coincides approximately with the harmonic one. 
To get rid of fast oscillations, one applies the method of averaging js^: One inserts for the 
amplitude x(t) = A(t) sin(6o) = A(t) sin(co'ot) and collects terms of the different time scales 
1 and ujq. So, one obtains an equation for the slowly varying amplitude A{t), and for the 
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phase 00 here for the van der Pol oscillator. 

A = 1/2 (^a - \(3A^^ A , (2) 
Qo = uJo . (3) 

Eq. ini describes the slow relaxation of the amplitude towards the steady state. A limit 
cycle exists for \A\'^ = The solution of the phase equation allows for the addition of an 
arbitrary constant phase. Thus, one can shift the state along the cycle without changing the 
energy. This allows for the adjustment of the phase, when the system is driven, or coupled 
to another oscillator. This adjustment is exactly what we understand as synchronization, 
or mode- locking. 

A Driven, Single Oscillator: If the system is driven externally by a harmonic force with 
angular frequency Ui ~ Uq, two time scales are present in the system, a fast one tj = ^ ~= 
and a slow one tg = ^^J^, ,^^ , and tg <^tf. The corresponding equation 



UJQ ' LOl—UJQ ' 



df 

X + lUqX — X— = uj-i^Rcos{uJit) , (4) 



We want to investigate now the dependence on the slow time and average over the fast one. 
To do so, we use the ansatz x(t) = A(t) sin(co'ii(: + 0), with the slow phase (It is a bit easier 
working with a complex amplitude, substituting x = Ae^'^^'^^'^, x = Auie^'^^^'^f}. By using 
again the averaging method one obtains after a few straightforward steps 

1 ^ 2\ UJlR 



2y^~4^^)^ ^cos0, (5) 

cP = ^ + ^-sm0. 6 

The second equation is known as the Adler equation Q,!^!; the general form for the phase 
equation depends on the nonlinear function f{x) and can be seen as the difference of the 
equations for the two fast variables ui, and uq: 



®' '''' = 6)2 - Gi = Ao^ + eg(0) (7) 



2 2 

with a 27r-periodic function q, and the two parameters detuning Aou = ~ u;o — cJi 

and locking term e = , which is proportional to the driving strength. To zero order. 
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one can assume a ~ constand the phase equation effectively decouples from the amplitude 
equation, which in turn is driven by the phase, similar holds for a first order (l9|. The 
Adler equation {q{4>) = sin(0)) has two stationary solutions, = 0, for |Ati;| < e. The 
stability of these fixed points, (ps and 0„, is determined by linear stability analysis. One 
puts (f) = (ps/u + ^'Pi inserts this expression into ((Zj) and Taylor-expands the expressions. 
The stable point is given by cos0s < 0, the unstable one by coscpu > 0, corresponding to 
exponential decay or growth of the perturbation 6(f). In general, if higher harmonics enter 
in sev Fig.|21 




FIG. 2: Stable (filled circle) and unstable (open circle) fixed point Qg, for the Adler equation 
to illustrate the mode locking. If detuning or coupling is varied, the phase is adjusted accordingly 
by the system to the stable fixed point. For Auj — e = both fixed points merge; this results 
in a saddle-node bifurcation, as sketched in the inset; the straight/dotted line corresponds to the 
stable/unstable fixed point. 



The transition to synchronization: If one varies the detuning from zero to e the phase 
varies from zero to (at most) tt. At the merging -, or bifurcation point, Aumax ~ e = 0, and 
the stable and unstable fixed point annihilate (opaque in Fig. |21). This type of bifurcation is 
known as saddle- node bifurcation js^. In the case of an unstable saddle one attracting and 
one repelling direction to/from the periodic limit cycle exist, for a stable cycle two attracting 
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directions exist, as for a stable node. 

Inside the synchronization region, the phase difference (j) = (pQ = const, and the phase 
of the observed oscillation is ujit + ^q. Let us turn to the phase dynamics outside the 
synchronization region. For 7^ 0, one can formally integrate Eq.Q to obtain 



t 



cmi/y ) — Alj 



(8) 



The period T = 27T/uJb results with the integration bound (p = (po + 2tt. The total phase 
rotates uniformly (f = uit + = uit + Ubt = AQt and one recognizes the typical beat 
phenomenon, as for the superposition of two harmonic waves; only that here the beating 
is strongly nonlinear. For harmonic oscillations the phase-slip of the beating is distributed 
over the whole interval. Close to the bifurcation from the synchronized region, one observes 
long epochs of nearly constant phase (p — <Pmax broken by relatively short intervals in which 
the phase rotates by 2tc. We observe this behavior as well for coupled organ pipes, this is 
illustrated below by a plot of the measured time signal close to the synchronized region. 

For a saddle-node bifurcation, one expects a square root dependence of Afl on the de- 
tuning Aa;. This indeed by expansion of the denominator in (|S)): 



-1 



|Afi| ~ 27r / — , (9) 

7-00 e/2q"{(pmax)(p'^ - (Au; - Auj^ax) 

= Ve|g"(0max)| ■ (Au; - AUmax) - a/ AcJ - AUmax (lO) 

A very detailed description, including many examples and generalizations is given in jiol . 

Two coupled oscillators: Now, we will turn our attention to the case of two oscillators 
and focus to the equations for the phases. From the above, it is now clear that any pair of 
uncoupled, self- sustained oscillators close to a limit cycle can be written in terms of phase, 
(61, 62), and amplitude, {Ai, A2), in the following form: 

= (11) 

A, = -7 (Ai - Afl) , (12) 

with i = 1,2. For weak coupling, one can again perform an expansion and apply the method 
of slowly varying amplitude such that the phase equations are, analogous to ^ 

6)1 = cui + eGi(ei,e2) , (13) 
6)2 = a;2 + 6^2(61, 62) . (14) 
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For uji ~ ulI2 the phase difference = 0i — 02 is a slow variable. The discussion can now 
follow the above, with the difference that now the oscillators have to adjust themselves 
mutually, whereas above one oscillator adjusted its frequency to an external driving. 

We have to distinguish very carefully and clearly between the single angular frequencies 
oji, measured for a decoupled system, and the angular frequencies Bj measured for the cou- 
pled system. For a clear nomenclature, we call the difference of the uncoupled frequencies 
detuning, and denote it by the symbol A/ = {2ti)~'^ {ijJ2 — uji) . The frequency difference oi the 
coupled system is always denoted by Az/. A typical plot visualizes Az/ against the detuning 
A/. One observes a plateau Az/ = for a synchronized region and a linear dependence for 
uncoupled oscillations. Below, we will use the symbol / for uncoupled frequencies, and u 
for measurements of the coupled pipes, analogously we use A/ and Az/. 

At the end of this section, we would like to hint to two facts: i) the theory presented 
above holds only for small deviations from harmonic oscillations and for small coupling. 
For large parameters, in principal the considerations still hold, but quantitatively deviations 
are expected, ii) Even though it seems appealing to explain the frequency locking of organ 
pipes in such a simple way, it is not completely satisfactory because organ pipes are extended 
systems and such a simple description does not take into account the aeroacoustics which is 
eventually needed for a complete understanding. 

III. MEASUREMENT AND RESULTS 
A. Measurement Setup and Signal Analysis 

Measurements were carried out on a miniature organ especially made by Alexander Schuke 
GmbH for physical investigations. Just like a real organ it consists of a blower connected 
by a mechanical regulating- valve to the wind-belt and finally a wind-chest upon which the 
pipes are positioned. Since only stationary behavior was of interest we joined the two pipes 
directly to the wind-belt with flexible tubes. The wind pressure was set by the organ builder 
to a mean value of 160 Pa. Over 10 s we measured a standard deviation value of 6 Pa in the 
pipe foot. The wind may be considered as very stable. 

The pipes used are stopped, a detailed description of the pipe-geometry is given in table 
H] The resonator top of the pipes was movable in order to adjust resonator- length and pitch. 
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To minimize losses, a gasket made of felt was applied. The pipes as delivered from the 
manufacturer were tuned to e and f in German notation, that is E3 and F3 in American 
notation. The applied driving pressure roughly yields a Reynolds number Re = 435. 

The acoustic signal was measured by a B&K 4191 condenser microphone positioned 
at the centerline between the sound sources. Further information on positioning of the 
acoustic sources and the microphone is given in sections IlII Bl and IIII CI Measurements took 
place inside the anechoic chamber of the Technical University of Berlin. The rooms size 
is W ■ L ■ D = 20m ■ 17m ■ 7m giving a volume of 2000 m'^. The low cut-off frequency of 
this room is 63 Hz. Temperature, humidity and pressure conditions were kept stable to keep 
there influence on the measurements negligible. 

The recording of the sensor signal as a function of time was achieved via hard-disk 
recording with a M- Audio 2496 Soundboard at a sampling rate of 44.1 kHz and a resolution 
of 16 bit and saved in PCM format. Data analysis was performed with Matlab. All 
Spectra shown in the results are Fourier transformed out of 10 s (441 000 samples) of time 



TABLE I: Pipe geometry. 



pipe-body 


overall length (without foot) 


530 mm 


clear length (adjustable) 370 . 


. . 450 mm 


wall thickness 


6 mm 


cross section 


shape 


quadratic 


clear width 


37 mm 


clear depth 


48 mm 


mouth 


cut-up height 


7 mm 


cut-up width 


37 mm 


flue-exit height 


0.4 mm 


foot 


length 


320 mm 


toe hole diameter 


15 mm 
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FIG. 3: Sketch of the experimental setup. The microphone is in plane with the pipe mouth, the 
view is at an angle from above. 



signal, giving a frequency-resolution of 0.1 Hz. During Fourier transformation a symmetric 
Hann window was applied to time data. All spectra have been averaged. From these spectra 
the amplitude of the harmonics was read out by a routine searching for the maximum signal 
level within a frequency range of ±0.3 Hz. 



B. Organ Pipe driven externally by a Loudspeaker 

To investigate the influences of a external sound source onto the organ-pipe we posi- 
tioned ax loudspeaker directly beside a single organ-pipe. The loudspeaker was a 2- Way 
Active System MS 16 from Behringer. The woofers cone-diameter is 10 cm. The speaker 
was connected to a HP 33120A waveform generator, delivering a sinusoidal signal. The 
sound pressure level emitted by the loudspeaker was set equivalent to the volume of the 
fundamental organ tone within 81.5 ± 0.25 dB. At this level the distortion of the speaker 
was at THD=1.5%. The experiment started with the pipe at a fundamental frequency of 
168.3 Hz and the loudspeaker at 165.8 Hz. During the measurement the organ-pipe remained 
unchanged but the frequency of the loudspeaker was very slowly, every 40 seconds, increased 
by 0.1 Hz. This gives the possibility to average ten seconds of signal four times to suppress 
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noise. The microphone was positioned at the centerhne between pipe and loudspeaker at a 
distance of 15 cm. A picture of the setup is given in Fig. |3 




FIG. 4: Fotograph of the measurement setup. Left: loudspeaker and pipe sounding together, right: 
two "coupled" organ pipes 

The synchronization can be observed clearly by the plateau in the plot of frequency dif- 
ference, Au against detuning. A/. As well the predicted square root function is observed 
at the transition to synchronization, cf. Fig. El The amplitude recorded with the micro- 
phone shows a positive and negative peak, at the beginning and end of the synchronization, 
respectively. Invoking the explanations of synchronization theory, cf. Eq. d this amounts 
to a change of the relative phase from to vr. If loudspeaker and pipe are considered as 
monopole sources, their field at the centerline, where the signal is recorded, can be obtained 
as a superposition of two fields with equal, measured amplitude and phase difference (po, 
decaying as 1/r: A{r) = Aq/v + A^/r cos (po (remember that we adjusted the sound pressure 
level of loudspeaker and pipe, so Agpeaker = Apipe = Aq). For (po = interference is positive, 
00 = TT implies negative interference. This is exactly seen in Fig. Eb, where we plot the 
peak amplitude of the loudspeakers frequency (open circle), and the first harmonic of the 
organ pipe (filled circles), respectively, over the detuning. For zero detuning, 0o = 7r/2. 
This corresponds to an interval of vr within the synchronization region; i.e. from positive 
to negative superposition. One notices that the phase (po = 7r/2 right at A = Aq. This 
corresponds to an equal distribution of energy at this point. The increase of the amplitude 
at the loudspeakers frequency at the left hand of the synchronization region comes from 
the above mentioned effect that the pipe sounds for long periods with the frequency of the 
speaker to undergo a fast phase slip, when close to be synchronized. This results in a higher 
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FIG. 5: Upper left (a): Frequency locking of organ pipe and loudspeaker. Top: observed frequency 
difference Av ascertained with pipe and loudspeaker sounding together versus detuning Af (open 
circles). In the synchronization region, a very clean plateau is observed. A saddle- node bifurcation 
(filled diamonds) occurs at the edges, as predicted by Eq. Q. Bottom left (b): Measured amplitude 
signal in dB. The phase varies from to vr over the synchronization plateau. The region is asym- 
metric; Av G (—0.3, 1.1) within the frequency resolution. The circles indicate the peak amplitude 
at the loudspeakers frequency, and the pipes frequency (squares). The superposition of the signals 
including the phase shift along the plateau is plotted by filled diamonds, the agreement is quite 
good. Right (c): absolute measured signals for pipe and loudspeaker sounding together {v, open 
circles) and the signals for speaker alone (squares) and pipe alone (straight line). The detuning 
is indicated by Af, as the difference of the "uncoupled oscillators", pipe and speaker alone; the 
frequency difference, Av results from the subtraction of the measured signal for pipe and speaker 
sounding simultaneously, coupled and the varied frequency, in our case of the loudspeaker. 



C. Two Pipes standing side by side 



To measure the dependence of frequency locking on the detuning of the two pipes, both 
pipes were installed on a common bar directly side by side (see Fig. HJ One of the pipes has 
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been kept at the fixed frequency of 170.1Hz, the other pipe was detuned in variable steps. 
With two pipes connected to the wind-chest one cannot detect the fundamental frequency 
of one pipe by simply turning of the other through closing the pipes valve. This would 
raise the wind pressure for the investigated pipe and consequently also the frequency itself. 
In order to acquire the value precisely the pipe at a fixed frequency was made soundless 
by putting a little absorber-wedge into the mouth. Doing so the air-sheet oscillation and 
therefore sound generation is suppressed, but the air consumption remains the same, keeping 
the wind-pressure stable. 

In H, the frequencies z/i, 1/2 of the coupled pipes were plotted against the frequency 
detuning of the uncoupled pipes. A/. As explained above, a clearer characterization is 
achieved when the frequency difference Az/ is plotted versus A/ [3, 1^- We plot this 
quantity in Fig. IHl Clearly, the typical behavior of frequency locking can be observed. As 
synchronization theory predicts, the bifurcation close to the ends of the locking region is of 
saddle node type (filled diamonds in Fig. IHt)- 

Now, stressing the analogy with two c oup led oscillators, a change of the relative phase 
from (j)Q — Tx/2 to 0o + 7r/2 is expected |l2| (0 = 27r/). This has consequences for the 
amplitude of the emitted sound. The radiated sound wave in the far field can be roughly 
approximated by a spherical wave, and the field at the measurement point r is obtained 
by superposition of the two waves emitted at the pipe mouths positions. Assuming the 
sources to be point-like and considering that the pipes are situated directly side-by-side, the 
amplitude along the center line is A = 2Ao/r (1 + cos0o), this has already been described 
in |l^. The notation is as above, please note that here Apipei = Apipe2 = Aq, because both 
pipes have approximately the same amplitude, cf. Fig. IHl 

Varying the frequency mismatch A/ we observe a gradual increase of the amplitude at 
the centerpoint, where the microphone is situated. In Fig.lHb, the result of the measurement 
is displayed: We see an excellent agreement of our estimate with the measurement. This 
behavior is observed as well in the spectral plot. Fig. [3 where one can observe the collapse 
of all sidebands in the measured spectrum to one single frequency as a very sharp transition. 

For higher harmonics, the phase difference by 0o does not imply destructive interference 
because they are slaved and follow the first harmonic. In other words, the relative phase 
for the nth harmonic in the synchronization region lies in the interval {n{(j)Q — TT/2),n{(j)Q + 
7r/2)). For the second harmonic, this implies in-phase oscillations at A/ = 0. This is heard 
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FIG. 6: Frequency locking of two organ pipes. The plot shows the observed frequency difference Aj^ 
versus the detuning of the uncoupled pipes, Af (circles). Upper left (a): In the synchronization 
region, a very clean plateau is observed. As in the case of external driving, the saddle-node 
bifurcation is observed at the edge to the synchronization region (filled diamonds). The agreement 
is very good. Lower left (b): Amplitudes at the first harmonics of the pipes (circles, variable- 
frequency pipe;squares, fixed-frequency pipe) against detuning; the sharp decrease of the amplitude 
at Af = indicates an antiphase oscillation at the pipe mouth. The analytical curve obtained 
using synchronization theory is shown with filled diamonds, assuming (po = vr — the agreement is 
excellent. Right (c):Plot of the absolute frequencies (circles/squares as above),, the straight line 
shows the 170.1 Hz of the fixed-frequency pipe. 

acoustically by a dominance of the octave around Af = 0. At the edge of the synchronization 
region, however, there is a phase difference of vr and destructive interference is expected. 
This is confirmed well by the quantitative measurement of the amplitude dependence of the 
second harmonic, shown in Fig. |S1 The agreement with the estimate from superposition 
of two monopoles is still good, but deviations occur. Higher harmonics are in qualitative 
accordance with the above ideas, but quantitatively more and more differences are observed. 
This can be expected, because the approximation by a monopole source does not hold. 

When two pipes are coupled diffusively, a so-called oscillation death can occur, leading 
to a complete silence of both pipes. The whole input power would then be converted to 
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heat without radiation. We thus observe a strong impact of the synchronization of the 
pipes on the acoustic properties of the emitted sound. From our results, we can rule out an 



oscillation death for two coupled pipes with zero detuning, as suggested in 



la]. Rather, the 



two oscillators radiate two antiphase waves; this results in a vanishing acoustic signal. 



D. The Coupling of the Pipes 

As in the case of the external driving by a loudspeaker, the two pipes interact by the pres- 
sure perturbation generated at the pipes mouths. But now, there are two sources of pressure 
contributions: the acoustic and the aerodynamic part from the vortices itself. Which one 
is dominant depends on the ratio of their amplitudes. In the case of the organ pipes, the 
acoustic pressure is very large due to the selected amplification of the resonator frequencies. 
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FIG. 7: A sharp transition to synchronization is observed and the sidebands from beating in the 
measured spectrmn collapse to a single frequency in the synchronization region. The amplitude is 
encoded according to the levels on the right. 



17 



However, the wind field at one pipe can literally "blow away" an adjacent vortex. On the 
basis of our observations, we can only state that definitely, the acoustical field is sufficient 
to give raise of synchronization. The influence of the wind field must be subject of careful 
future work including flow visualization on the scale of several vortex diameters. At the edge 
of the synchronization region, there must be a phase difference of 0o =t '?i"/2, which corre- 
sponds to a shift of a quarter of a jet oscillation. The jets then undergo a mutual reordering 
into a complex spatio-temporal three dimensional pattern. Again, the observation of the 
pattern was beyond the facilities of the current experimental setup. 

All our observations are consistent with an interpretation in the frame of synchronization 
theory of two oscillators. A full description of the physics has to take into account the 
origin of the sound generation — the oscillations of resonator and the vortices generated by 
two jets. Most likely, there is a competition of the interaction between jets and resonators. 
Which mechanism wins depends then on the ratio of the respective amplitudes. A comment 
should be made on why the synchronization description by two oscillators holds so well. In 
the measured Reynolds number regime, there might be three-dimensional patterns at the 
pipe mouth, if the aspect ratio is large enough. To check this in our particular case requires 




FIG. 8: Amplitude and frequency difference for the second harmonic. Here, the phase is doubled 
and an annihilation occurs right at the edge of the synchronization region, seen as a sharp decrease. 
The curve predicted by theory shows deviations for the second harmonic. 
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detailed experiments. Neglecting three dimensional structures, and assuming that the air 
sheet is homogeneously oscillating, a description by oscillator models like Stuart-Landau or 
a van-der-Pol equation could be sufficient to describe the two dimensional oscillations and 

n 

the transition to turbulence Since the main sound production is due to the interaction 
with the labium, an oscillator model suits well. From synchronization theory it follows 



that any self-sustained oscillators will generica^ 



coupled, leading to the observed behavior 



ly follow the synchronization scenario when 



lol 37I . The simulations in are completely 



in agreement with this fact. We are convinced that full understanding may be achieved by 
studying the details of the coupling mechanism. 

To investigate the sensitivity to coupling strength, we measured the frequency spectrum 
about the ffist harmonic while varying inter-pipe distance (measured between the outer 
walls of the pipes). The detuning was fixed to 0.7 Hz. The corresponding plot is depicted in 
Fig. El For large distances, the typical interference pattern of two noninteracting oscillators 




FIG. 9: Frequency spectrum of the coupled pipes in dependence on the distance, given in the 
legend. The pipes decouple more and more when farther apart. For small distances, a sharp peak 
at the synchronization frequency is observed, for large distances, the peak broadens more and 
more, at very large distances the spectrum of the uncoupled pipes is recovered with the typical 
beat phenomenon. To distinguish different graphs, an offset has been added to each curve. 
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is observed. As the pipes come closer to each other, the individual, sharp peaks of half width 
^ 0.1 Hz reduce in amplitude and at the same time the spectrum broadens to a "hump" 
with half width ^ 1 Hz. Peaks typical of the beating from linear superposition sit on the 
hump. Coming even closer, full synchronization is observed with one single, very sharp peak 
of again half width ~ 0.1 Hz. That means the hump is about ten times broader than either 
of the peaks for the uncoupled system or the synchronized one. 

The basic observations can be explained by synchronization theory. In Eq. El the square 
root behavior of the frequency difference at the transition to synchronization (the bifurcation 
point) is derived. Lowering the coupling strength is equivalent to shift the bifurcation point 
towards zero. At a certain distance the detuning leaves the synchronization region and 
one should observe a (strongly nonlinear) beat. The peaks of the beat are observable as 
small side peaks on top of the broad peak. For a pure beat, only sharp peaks should be 
observed. Here, the duration of the large time interval -we do not want to call it period- 
varies slightly such that the sidebands wiggle a bit, altogether this leads to a broad peak. To 
illustrate this behavior, we plot in Fig.llOlthree different time signals for the distances 10 mm 
(synchronized), 15mm (shortly before synchronization), and 25mm. The fast oscillations 
look very similar for all of the distances, as exemplary given by Fig. ^1 (bottom, right). 




Tim=(.) TimB(m.) 



FIG. 10: Time signal for three different distances. Upper left: 10mm, lower left: 15 mm, upper 
right:25 mm. The lower right graph shows a typical fast oscillation as found for all three distances. 
The transition from synchronization to beat happens between 15 and 25 mm. For 10 mm the am- 
plitude is constant and the phase rotates with the fast frequency. At 25 mm the beat phenomenon 
is observed with a constant amplitude over relatively long times (e.g. ) and a quick phase slip (). 
At 15 mm, the system shows some intermittent bursts at irregular times. 
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A complete understanding why such a broad hump appears could not be achieved. There 
are several possible explanations. According to classical synchronization theory the fre- 
quency fluctuation cannot be explained but by noise, e.g. in the air supply in the wind 
system. Another option are the small fluctuations in the vortices at the pipe mouth which 
can be sufficient to cause the very small variation in the beat frequency. Finally, one can con- 
sider the bifurcation to synchronization as a kind of phase transition, where typically large 
fluctuations occur, including chaotic behavior. By means of our measurement we cannot 
decide which of these scenarios is true. Finally, we want to mention recent results concern- 
ing whistling of two nearby Helmholtz resonators show a similar transition in the phase of 
the acoustical signals emitted by the resonators 0]. Interestingly, the phase relation is lost 
at distances of the order of the ones we find (25mm, if the wall thickness of the pipes is 
included). Thus, this observation might be explained by synchronization as well. 



IV. CONCLUSION 



We presented measurements on an external driving of an organ pipe by a loudspeaker and 
on the mutual influence of two organ pipes. The observed behavior is completely consistent 
with an explanation in the frame of synchronization theory. For the dependence of the 
synchronization on the "coupling strength", measurements determining the dependence of 
the frequencies on the inter-pipe distance have been carried out. They reveal a broadening 
of the peak at the synchronization frequency with increasing distance. This broadening can 
be explained by quick phase slips as predicted by theory close to a saddle-node bifurcation. 
The smearing out of the side peaks of this beat can be either due to inconstant wind supply 
or due to the fluctuations caused by the oscillating air sheet at the pipe mouth. To clear 
this question further measurements have to be carried out. 

From an acoustic point of view, one can address the question of how to position two 



organ pipes close in frequency, 
and error in the last centuries 



his has been intuitively solved by organ builders by trial 
17| . Our work might give quantitative hints on how large the 
inter-pipe distance needs to be to suppress mutual influence, and on details of the coupling 
mechanism. For example, avoiding an amplitude minimum for the first harmonic is highly 
desirable. 

From an aerodynamical point of view, the above scenario requires more detailed inves- 
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tigations to understand the full dynamics of the coupled resonator/jet system, although a 
model consisting of two mutually coupled oscillators seems to be sufficient for all qualitative 
questions. In addition the more involved setup of more than two organ pipes is an interest- 
ing subject for further investigations. From our measurements, one can clearly say that the 
pipes do not show an oscillation death; rather, antiphase sound radiation yields the observed 
weakening of the amplitude. 
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